Detection of Pair-Superfluidity for bosonic mixtures in optical lattices 
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We consider a mixture of two bosonic species with tunable interspecies interaction in a periodic 
potential and discuss the advantages of low filling factors on the detection of the pair-superfluid 
phase. We show how the emergence of such a phase can be put dramatically into evidence by 
looking at the interference pictures and density correlations after expansion and by changing the 
interspecies interaction from attractive to repulsive. 
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Ultra-cold atoms in optical lattices are presently one 
of the best environments for the study of exotic quan- 
tum phases 0, The experimental demonstration of 
the superfluid to Mott transition has opened the way 
to the study of strongly correlated phases in lattices. In 
the case of mixtures of different atomic species or differ- 
ent internal levels, new phenomena related to quantum 
magnetisms and spin physics arise [3]. 

The quantum phase that is central to this work is the 
pair-superfluid (PSF) phase for a mixture of two different 
bosonic species. This phase consists in the formation of a 
superfluid of pairs where atoms of different species prefer- 
entially hop together in the lattice. It is characterised by 
the vanishing of the single-species order parameter and 
by the emergence of a pair-order parameter. This prob- 
lem is related to the historical question formulated, e.g., 
by Nozieres and Saint James [5( , concerning the impossi- 
bility of coexistence of single and pair superfluidity. This 
impossibility lies in the fact that a single species con- 
densate is associated with the macroscopic occupation 
of the zero momentum (quasi-momentum) state for each 
species, while a pair superfluid requires the macroscopic 
occupation of the zero relative momentum, without any 
constraint on the value that the single species momenta 
can separately take. As a consequence, single species 
condensation is destroyed. In this work, we will see how 
this picture is directly reflected into physical quantities 
accessible in experiments. 

The existence of a pair-superfluid phase in lattices has 
been already pointed out in several papers [1, 0, H, H, OH 
HP . HH, [IH . The main emphasis has been put on the situ- 
ation of equal density for the two species leading to total 
integer and half integer filling factor. For integer filling 
factors, PSF arises in the regime where the interspecies 
interaction almost completely compensates the repulsive 
intraspecies interaction. Unfortunately, the precision on 
the values of the interaction strengths and the very small 
values of the tunneling parameter required to get pair 
superfluidity, makes this phase almost unaccessible ex- 
perimentally. At half filling factor for each species (total 
integer filling), the PSF phase (analogous to the x — y 
ferromagnet) is predicted in a larger region of the phase 
space. At low tunneling and in the presence of asym- 
metries between the interaction and tunneling parame- 
ters of the two species, the PSF phase competes with 



the insulating-like anti-ferromagnetic ordering. Instead 
for total incommensurate filling factor, no insulating-like 
phases (Mott or antiferromagnetic-like) exist. In this 
regime, easily accessible signatures for the experimental 
observation of pair superfluidity are available. In this pa- 
per, we would like to complement the predictions in [l3(, 
discussing the role played by the two-body momentum 
distribution and commenting on the effect of interactions 
in the expansion. 

At total incommensurate filling factor and zero tem- 
perature, two important phases are naturally conceived 
[14| : (i) a double superfluid (2SF), where both species 
are independently superfluid and single-species coherence 
exists, and (ii) a PSF phase, characterized by pair co- 
herence. Assuming that interactions do not affect sig- 
nificantly the expansion of the atomic cloud after re- 
lease, all required information needed to distinguish be- 
tween the two phases are included in the pictures of the 
two species after expansion: first of all, the interference 
fringes, typical of single-species coherence, will appear for 
the 2SF phase and vanish in the case of PSF. Moreover, 
the density-density correlations between the two species 
after expansion carry information about the correlations 
in momentum space before expansion, which are dramat- 
ically different for the two phases. 

We consider two bosonic atomic species in a lattice, 
described by the Bose-Hubbard Hamiltonian 
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where a = a,b indicates the two bosonic species, 
a,-, hi, at, b\, and nf are respectively the annihilation, cre- 
ation operators and the density of species a and b at site 
i. The notation (ij) represents nearest neighbors. The 
intraspecies on-site interactions U a and tunneling param- 
eters J a depends in the standard way on the optical lat- 
tice potential and scattering lengths. Generally speaking 
the most favorable conditions for PSF are given by a com- 
plete symmetry between the two species, as far as interac- 
tion, hopping and density are concerned. This is the sit- 
uation that we will assume in this paper (N a = = N, 
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J a = Jb = J, U a = Ub = U), focusing on the possi- 
bility of changing the interspecies interaction U a b from 
negative to positive by tuning the interspecies scattering 
length via a Feshbach resonance over a wide range, as 
demonstrated in (l5j . 

The solution of Hamiltonian JT]) is a non-trivial many- 
body task. Quantum Monte Carlo (QMC) calculations 
provide the most accurate results 0, @, E3] • Evidences 
of the PSF phase have been recently obtained also by 
matrix-product-state 0, Il2l . [l3j and dynamical mean- 
field approaches [ll|. Instead the standard mean-field 
approximation based on a generalized Gutzwiller Ansatz 

*) = IL J2 n ,m fn,m\n a , rib)i, which in principle includes 
correlations between the two atomic species, misses the 
pair superfluid phase due to the impossibility of correctly 
accounting for second order hopping which is at the ori- 
gin of the PSF phase. Recently, a treatment based on 
a mean-field analysis of the effective Hamiltonian in the 
pair-subspace applied to a bilayer system of 2D dipolar- 
lattice bosons has proven successful to describe the PSF 
and pair-supersolid (PSS) phases [l~7l ]. 

In order to capture the basic physics of the pair cor- 
relations underlying the emergence of the PSF, in this 
work we employ a toy-model based on the exact diago- 
nalisation of ([T]) for a system of few atoms occupying few 
lattice wells in ID with periodic boundary conditions. Of 
course sharp phase transitions are not accounted for by 
our treatment, but we believe that the main conclusions 
remain valid also for larger systems and higher dimen- 
sions. 

The comparison between filling factor equal to and less 
than 1/2 is useful, since for filling factor exactly equal to 
1/2, there exist a particle-hole symmetry between posi- 
tive and negative interspecies interaction in the almost 
hard-core limit \U a b\ *C U, leading to the pairsuperfmid 
phase (PSF) for attractive interactions (pairing of two 
atoms of different species) and the so-called countersu- 
perfluid phase (SCF) for repulsive interactions (pairing 
of one atom with a hole of the other species). Instead, for 
equal filling factor less than 1/2 for each species, PSF still 
persist, but SCF pairing is suppressed. This different be- 
haviour for positive and negative interspecies interaction 
might help identifying the formation of PSF as discussed 
below. 

In order to put into evidence the differences between 
filling factor v — 1/2 and v < 1/2, we will consider a 
system of 4 wells and 2 atoms of each species (N w = 4, 
N = 2) and a system of 6 wells and 2 atoms of each 
species (N w — 6, N — 2), which are the smallest cases 
including the possibility of having equal filling factors 
v < 1/2 and non trivial on-site interactions. We look at 
the following quantities fl6| : 

V 2S f = (otoi+x) = (blb i+1 ), (2) 
Vpsf = (atfeja i+1 6 i+1 ) - (a\a i+1 ) (b\b i+ i) , (3) 
Vscf = {a\b iai+1 b\ +1 ) - (a\a i+1 ){bib\ +1 ). (4) 
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FIG. 1: (Color online) Phase diagram for N w = 4, 6 (upper 
and lower row, respectively) and N = 2. (a,d) Single-particle 
coherence N W V2SF, as denned in Eq.([2]); (b,e) Pair coherence 
N w Vpsf, as defined in Eq. ([3J ; (c,f) Counter-pair coherence 
N w Vscf, as defined in Eq.(j4|. The factor N w allows a better 
comparison between the two different lattice sizes. The re- 
gion of parameters of strong attractive interaction U a b < — U 
corresponds to collapse in a large system. 

The quantities V2SF, Vpsf and Vcsf characterise re- 
spectively the 2SF, PSF and SCF phases. We also con- 
sider single particle and two-particle momentum distri- 
butions in order to make a useful link to the experiments. 

In FigJTJ we show the phase diagram as a function of 
J and U a b- One can clearly identify the PSF (or SCF) 
regions as the dark regions in FigQ](a,d) where single- 
species coherence vanishes, and at the same time pair- 
coherence (or counter-pair coherence) is different from 
zero, namely the light regions in FigQJb,e) ( or FiglU c ) 
for SCF). As explained above, for v = 1/2, PSF and 
SCF are found respectively for attractive and repulsive 
interaction. Instead, for v < 1/2, SCF is absent. Fol- 
low filling, the crossover from 2SF to PSF is governed 
by the competition between the single particle hopping 
J and the energy cost for breaking a pair, equal to |/7 Q b|. 
For that reason, PSF is found for attractive interspecies 
interactions at sufficiently low tunneling parameter J ^ 
\U ab \. 

An important quantity accessible in experiments, 
which would provide an unquestionable proof of PSF, 
is the measure of correlations in the momentum distribu- 
tion, reflecting the fact that in the PSF and SCF phases, 
two atoms of different species form a pair and condense 
in the state of total quasi-momentum q a ± qb = (re- 
spectively for atom-atom and atom- hole pairs), as shown 
in Figj2ja,c,f). In the case of two non interacting su- 
perfluids (U a b — 0), the two species have completely un- 
corrected momentum distributions, i.e. n Ka ' b ' > (q a ,qb) — 
n^ a \q a ) x n^(qb), where n( a \q a ) and n^(q b ) separately 
present interference peaks at even multiples of the Bragg 
momentum qs (see Figl2jb,e)) i _as happens for standard 
single component condensates [3, EH • In the presence of 
interspecies interactions U a b 7^ 0, correlations build up 
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FIG. 2: (Color online) Two-body momentum distribution 
n {a ' b) {q a ,q b ) for TV = 2 and N w = 4,6, i.e. filling factor 
1/2 and 1/3 (upper and lower row, respectively). (a,d) re- 
pulsive interparticle interaction U a b = 0.5(7; (b,e) vanishing 
interparticle interaction U a b = 0; (c,f) attractive interparticle 
interaction U a b = —0.5(7. In all pictures J = 0.01(7. Two 
Brillouin zones are shown for clarity. 
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FIG. 3: Single species momentum distribution n' CT '(q CT ) for 
N = 2 and N w = 4,6, i.e. filling factor 1/2 and 1/3 (up- 
per and lower row, respectively). (a,d) repulsive interparti- 
cle interaction U a b = 0.5(7; (b,e) vanishing interparticle in- 
teraction Uab = 0; (c,f) attractive interparticle interaction 
Uab = —0.5(7. In all pictures J = 0.01(7. Two Brillouin zones 
are shown for clarity. 



in a very different way depending on the filling factor 
and on whether the interactions are repulsive or attrac- 
tive. For filling factors exactly equal to 1/2, the situation 
is almost symmetric for positive and negative U a b upon 
particle-hole duality for the different species. The mo- 
mentum correlations are opposite in the two cases, as 
shown in Fig[^a,c), showing SCF and PSF respectively. 
The 2SF phase is recovered in the vicinity of vanishing in- 
terspecies interactions (Figj2]Jb)). For equal filling factor- 
smaller than 1/2, a 2SF is obtained both for vanishing 
and repulsive interactions (Fig|2{d,e)), since the filling 
factor of one species does not match the filling factor of 
the holes of the other. Hence, only negligible momentum 
correlations exist for U a b > 0. On the contrary, attrac- 
tive interactions lead to PSF and very strongly correlate 
the two different species (U a b/U = -0.5, Figf^f)). 

The correlations in the two-body momentum distribu- 
tion are strictly related to single-particle coherence and 
strongly affect the visibility of the single particle momen- 
tum distribution, as shown in Fig. [3) The presence of 
momentum correlations in n^ a - b \q a , qb) lead to a reduced 
contrast in n^(q a ). Hence, some signatures of the for- 
mations of the PSF/SCF phases are provided already by 
the interference in the single-particle expansion pictures, 
which is the most easily accessible experimental method 
[l3|]. For instance, at low enough tunneling, for v < 1/2 
interference is expected at U a b > 0, while it disappears 
(under exactly the same conditions) by tuning U a b to 
negative values. 

Most important, the single species expansion pictures 
carry information about the momentum correlations. In 
fact, as demonstrated in (l9| . the direct measurement of 
the momentum correlations can be performed by look- 
ing at the noise in the single-species expansion pictures 
[lOjIHi- Assuming first that interactions do not affect the 



expansion, the single-species densities after time of flight 
are given by n^ OF (r a = q a t/m a ) = n^(q a )- Hence, 
in the case of PSF, where the correlations are of the 
type q a + qb = 0, we expect the two expansion pictures 
to be correlated at correctly rescaled opposite positions 
and possibly at corresponding points in different Bril- 
louin zones. 

In usual experiments, interatomic interactions are not 
turned off during the expansion, and their effect is con- 
sidered to be negligible due to the higher energy scales 
involved in the problem. In the present case, it would be 
a safe procedure to tune at least U a b to zero just be- 
fore releasing the cloud. However, this might not be 
so easily achievable in an experiment. For this reason, 
we have performed simple numerical estimations of the 
effect of the interspecies interactions on the expansion 
for two atoms released from a 4- well lattice (see FigQ|. 
While the effect of relatively weak attractive interactions 
(\U a b\/J ~ 20, see Fig|I£a,b)), the two-body momentum 
distribution is hardly modified during the expansion, we 
have seen that the two-body momentum distribution can 
be affected by attractive interspecies interaction U a b, f° r 
interaction strengths leading to pairing (\U a b\/J ~ 60, 
see FigQJc,d)). This effect tends to create also some cor- 
relations along the diagonal q a — qb, but does not destroy 
the correlations typical of the PSF phase (at q a = —qb)- 
Hence, recovering the two-body density after expansion 
via noise correlation measurements, the distinction be- 
tween PSF and 2SF phases still remains very clear. More 
dramatic effects of the interactions during the expansion 
take place only for values of \U a b\/J which are much be- 
yond the estimated onset of the PSF phase transition. 

Small asymmetries in the Hamiltonian parameters of 
the two species do not destroy the PSF phase. Instead, an 
unbalance in the densities of atoms a and b hinders the 
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FIG. 4: (Color online) Effect of attractive interspecies inter- 
action Uab on the expansion. (a,c) Two-body density after 
expansion after a time of flight _E r tTOF ~ 12, where E r is the 
recoil energy, in the presence of interspecies interaction U a b', 
(b,d) Free expansion at the same time of flight (shown for 
comparison). The strengths of interaction are \U a b\/J ~ 20 
in (a,b) and \U a b\/J ~ 60 in (c,d) respectively. 



formation of the pairs. Unfortunately, we are not able 
to quantify the realistic effect of the unbalance due to 
the small number of atoms considered. However, one can 
think of experimental procedures to create a sample with 
almost exactly equal populations of the two species. For 
instance, one could start with a Mott insulator at unit 
filling for both species, and then tailor the optical poten- 
tials, introducing a second laser at half the wavelength, 
such to split each lattice well into two equal ones. Apply- 
ing this procedure along two dimensions would lead to a 
filling factor of exactly 1/4 for each species. This proce- 
dure should remain valid also in the presence of an exter- 
nal harmonic confinement, especially in the case of very 



low tunneling, where the central Mott region at unitary 
filling is dominant with respect to the outer superfluid 
shell. Alternatively, one can create a unitary filled Mott 
region for both species in the trap center and then release 
the harmonic trap till the desired filling factor is reached. 
This method would favor the creation of the pairs in the 
Mott phase, which can then become superfluid once the 
filling factor is made incommensurate. 

The physical ingredient on which pair superfluidity re- 
lies is the second order hopping of two atoms of the dif- 
ferent species at once. This is closely related to the ex- 
change interaction, whose observability has been recently 
demonstrated in [22[ . The fact that pair hopping is a sec- 
ond order process in J, where J is assumed to be small, 
might seem to be discouraging for the experimental ob- 
servation of the PSF phase. However, exact QMC sim- 
ulations of this problem [Icj, in line with the results of 
our toy-model, predict the transition between 2SF and 
PSF, at half integer filling and symmetry between the 
two species, to happen at J k 0.1 \U a b\- Exploiting the 
possibility of having U a b of the same order of U, this leads 
to relatively large values for the critical tunneling. On 
the other hand, careful analysis about the critical tem- 
perature and entropy for the formation of the PSF phase, 
as recently done in [23[, are required. 

Our model provides an oversimplified description of the 
system. We believe however that it includes correctly the 
fundamental ingredients of the physics involved. A more 
quantitative analysis based on exact numerical calcula- 
tions, including the effects of two-species unbalance on 
the formation of the paired phases, will be subject of 
future work. 
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